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1. Introduction 

The potential between two static colour charges in QCD is a quantity rich in features. At 
large distances r between the two charges the phenomenon of string breaking is observed for the 
theory with dynamical quarks [[[]]. The quenched theory is expected to asymptotically coincide with 
bosonic effective string theory for r — > oo [FJ]. At intermediate distances a scale ro ^ can be defined 
which is readily computed in numerical simulations and can be related to phenomenological models 
of quarkonium. This scale can then be used to set the overall scale of a simulation and/or set the 
relative scale for simulations at different lattice spacings. Finally, at short distances the shape of the 
potential is described by renormalised perturbation theory. In fact, by taking the second derivative 
with respect to the distance a renormalised coupling can be defined. In this way the static potential 
connects the non-perturbative regime and the perturbative regime and is also a quantity where one 
expects large effects of dynamical quarks, i.e., large differences between the quenched and the 
unquenched result. 

Here we report on an ongoing effort to measure the static potential on the configuration en- 
sembles generated by CLS (Coordinated Lattice Simulations) 1 . They were generated with the 
DD-HMC software package 2 , which implements two degenerated flavours of improved Wilson 
fermions and the Wilson gauge action. There are ensembles at three different values of j8 (i.e., 
three lattice spacings) and several values of the sea quark mass. The first objective is to provide the 
scale ro. It can be compared to other methods of scale determination pj|, used to compare dimen- 
sionless quantities among different collaborations, perform scaling analysis and preliminarily set 
the overall scale in physical units 3 . Since the data presented here are not based on all the available 
statistics and some details of the analysis might change, the results presented here are preliminary. 

The report is organised as follows. In the next section the techniques used to extract the static 
potential with low statistical fluctuations and small systematic errors are summarised (see [j|] for 
more details). In section || the scale ro is determined and extrapolated to vanishing sea quark mass. 
In section || applications beyond scale setting are collected, i.e., a comparison of the quark mass 
dependence of ro with other collaborations and the determination of a renormalised coupling. 

2. The static potential 

To determine the static potential V(r) we measure rectangular Wilson loops. They have exten- 
sion t and r in temporal and spatial direction respectively. Schematically, we measure 

C(t) = (tht), r fixed, (2.1) 

r 1 

where the brackets denote the expectation value with respect to the two flavour QCD measure 
and the square represents the product of link variables around a closed rectangular path. This is 
equivalent to a static quark-anti-quark pair at spatial separation r propagating distance t in time. 
For large euclidean time separations the signal is dominated by the ground state, which coincides 
with the static potential 

C(O f =°°A(r)e- y M', rfixed. (2.2) 

'https://twiki.cern.ch/twiki/bin/view/CLSAVebHome 
2 http://luscher.web.cern.ch/luscher/DD-HMC 

3 Because of the uncertainty coming from the phenomenological models used to give ro in physical units, this is 
somewhat unsatisfying and should be seen as an intermediate step. 



2 



Scale ro and the static potential from the CLS lattices 



Bjorn Leder 



fit of cff. mass 
o projection method ■ 
O GEVP, t a = 2 ■ 



_ 0.3535 - 

S 0.353 - 
I 

* 0.3525 i 



£ = 5.5 
= 5.3 



(/a 



Figure 1: Effective "mass" plot for the Wilson 
loop correlator matrix for r = 6a and to = 2a. 
See text for a detailed explanation. 



Figure 2: Static potential at two values of the 
lattice spacing. The static energy E s tat [§] is sub- 
tracted to obtain a renormalised quantity. 



A straightforward determination of V (r) from this definition with the original link variables of the 
configurations is bound to suffer from large statistical and systematic errors for two reasons. First, 
the signal to noise ratio will decrease exponentially for large r and t due to ultraviolet fluctuations. 



Second, the overlap with the ground state (encoded in A(r) in eq. (2.2)) will be poor. 

To remove the ultraviolet fluctuations the original link variables are replaced by smeared links. 
To this end one level of hypercubic smearing ^ with parameters a\ = 1.0, OC2 = 1.0, 0:3 = 0.5, 
referred to as HYP2 [f7|] was found to give the best result. The whole analyses was carried out 
with a second parameter set ot\ = 0.75, OL2 = 0.6, 0C3 = 0.3, referred to as HYP1. When not 
specified, the parameter set HYP2 is used. The smearing of the temporal links can be understood 
in terms of choosing an action for the static quarks and the smearing of the spatial links in terms 
of a redefinition of the operator that creates the static quark-anti-quark pair (see [||] for a detailed 
derivation). Using smeared links greatly improves the signal to noise ratio at large time separations. 

Since the exact wave function of the desired ground state is unknown the usage of a variational 
method is mandatory to improve the overlap. As already explained further smearing the spatial 



links on the left hand side of eq. (2. 1 ) is equivalent to a redefinition of the operator that creates the 
static quark-anti-quark pair. Intuitively, by smearing the generated sate becomes more and more 
extended. This way one obtains a correlator matrix 



C bn (t) = (tOt), 



r fixed . 



(2.3) 



In particular the indices l,m = 1, . . . ,M specify the n^ m levels of spatial HYP smearing with pa- 
rameters GC2 = 0.6, as = 0.3 4 that are applied to the spatial links. Throughout this report we use 
M = 4and m = l-\. 



From the correlator matrices (2.3) effective "masses" are extracted with two different methods. 



The starting point in both cases is the generalised eigenvalue problem (GEVP) 

C(t)\pa(t,t Q ) = X a {t,t Q )C(t Q )y a (t,t Q ), a = 0,...,M-l. 



(2.4) 



Only spatial links are involved, thus only two parameters are needed. 
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Figure 3: Chiral extrapolation of ro/a at 
j8 = 5.2. QCDSF data from [jl|] for comparison. 
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Figure 4: Chiral extrapolation of ro/a at 
j8 = 5.3. QCDSF data from Q for comparison. 



From the generalised eigenvalues one directly obtains 

E a {t + §,«)) = ln(Aa + *,*))) = + j 8„e- (£M - £ « )( ' + f ) + . . . , (2.5) 

where £o = V (r) and the higher states contributions are expected to die out with large Em — E a or 
faster. In [Q] the authors were able to prove this for t <2tQ. In fig. [I] Eo(t + f^o) 1S plotted for 
an intermediate r and to = 2a (blue circles) 5 . In order to quantify the contribution of higher states 



we perform a fit with the two terms of the right hand side of eq. (2.5) including data points for 
a = 0, 1 and t < 2?o (red dotted curve). The term modelling the higher states is used to estimate the 
systematic error due to them when extracting the ground state at a given time t (see below). 

Finally, the projection method of flTij ] is used to obtain a second determination of E a (t + 
|,fo). There the correlator matrices are projected to the ground state generalised eigenvector 
v = y/o(^o + a,to) 

f(t) = v T C(t)v. (2.6) 

The resulting numbers are locally fitted by a single exponential, i.e. three successive values f(t p ), 
t p = t-a,t,t + aare fitted to f(t p ) = beT^^v (black circles in fig. [I]). 

With the parameters as in this report we found the projection method to be slightly more stable, 
i.e. to exhibit longer plateaus. Therefore we extract the static potential from the black circle in fig. 
[I] with the smallest sum of statistical and systematic error (black dashed-dotted curve and grey error 
band). 

Applying this procedure to all values of r we obtain the static potential as a function of r. In 
fig. H we plot the result (made dimensionless by multiplying with ro, see below) for two values of 
the lattice spacing at roughly the same sea quarks mass. The plot shows that lattice artefacts are 
small. 



5 Statistical errors are determined using the method and program of [|lo|]. 
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Figure 5: Chiral extrapolation of ro/a at 
j8=5.5. 



Figure 6: Renormalised quantity c(r). Nf = 2 
data are at roMps « 1. See text for details. 



3. Scale ro 

The scale ro, introduced in is defined in terms of the static force F(r) = V'(r) by solving 

^(r)| r=ro = 1.65. (3.1) 



Its physical value is ro « 0.5 fm, thus it is sensitive to the non-perturbative character of the theory. 
Off the lattice it can only be determined through phenomenological potential models. On the lattice 
the static force is computed from the potential as the finite difference 



F(r I ) = \\V{r)-V(r-a)], 



(3.2) 



where rj = r — a/2 + 0(a 2 ) is chosen such that in ([3j2f) at tree level all lattice artefacts cancel out. 
To solve (3T) the force is then locally parametrised by a 2 F(r) = /o i^^/r 2 . Taking the two 



values that enclose the solution, the parametrisation is uniquely determined as well as r / a. This 
lattice definition of ro/a ensures that it is a smooth function of the bare parameters j3 and mo. To 
estimate the systematic uncertainties due to the parametrisation the same procedure is repeated 
with an additional term faa /r and three successive r-values. They are found to be negligible. 

We analysed three run sets at three different lattice spacings with /3 = 5.2,5.3,5.5. At each 
lattice spacing there are runs at 3-6 quark masses. For a scale setting we decided to extrapolate to 
the chiral point. In figs. 0-S the measured values of ro/a are plotted versus the subtracted quark 



mass am q = amo — am cr . Here we give the result for both parameter sets HYPl and HYP2. The 
extrapolation to am q = is done via a linear fit. The result of the fits are given in the legends of the 
plots. The HYP2 results are compatible with the HYPl ones and have the smaller error. Thus we 
list here for ro/a(fi) in the chiral limit: 

r /a(5.2) = 6.05(5) , r /a(5.3) = 7.05(3) , r /a(5.5) = 9.59(16) . (3.3) 

Using our new determination of ro /a we get a preliminary update on the A-parameter of [[7|] : 



r A 



Nf=2 
MS 



0.73(3)(5), 



(3.4) 
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Figure 7: Plot of dimensionless and physical quantities comparing different efforts to simulate two flavour 
QCD. The right panel is a blow-up of the region x = . . .2. 

where the first error comes from r^/a and the second from the running of the coupling. 

In figs. H and |] we also plot data from the QCDSF collaboration [|12|]. Since they are using 
exactly the same action, the results at coinciding bare parameters should agree within errors. The 
most tension is observed between their lightest point and the heaviest of this report at j3 = 5.2 in 
fig. ||. Also in fig. ^ the difference is hardly explainable by the difference in the /3-values (5.3 
versus 5.29). The most plausible explanation is that they used a global fit including many terms 
to model the potential (and thus the force) 6 over a wide range of r values. This way the statistical 
error in r^/a is reduced, but at the price of large systematic uncertainties. 



4. Physics results 

First, we present a quantity for which large effects of dynamical quarks can be observed. A 
renormalised, physical quantity can be defined in terms of the derivative of the static force 

c{r) = \r i F ! {r). (4.1) 

In [0J it was determined with high precision in pure gauge theory. With the data presented here we 
are able to compare the pure gauge case to the Nf = 2 theory. On the lattice we write c in terms of 
a finite difference 

c(r) = ^r 3 [V(r + a)+V(r-a)-2V(r)} , (4.2) 

where f = r + 0(a 2 ) is chosen such that at tree level all lattice artefacts cancel out. In fig. |6| we 
plot c for the three values of j3 at roughly the same sea quark mass (r$Mps ~ 1). For comparison 
we also plot the Nf = data [0]. The analytic curves in the plot are the 3-loop perturbative curves 
(dashed-dotted for Nf = and solid line for Nf = 2, spread due to uncertainty in the A-parameter) 7 , 
the universal value —n/12 from bosonic string theory JT3| , |l4] ] (that the Nf = data is approaching 
asymptotically for r — > oo), the value of c in the Cornell JTBj ] potential and the curve derived from 
the Richardson JT6[ ] potential. 

6 private communication with P. Rakow 

7 For details on the perturbative expressions for c(r) and how it is related to a renormalised coupling we refer the 
reader to 
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Second, we present in fig. [7] a comparison to other efforts simulating QCD with two light 
quarks. For the x-axis we define a dimensionless quantity x = r$Mp S (where Mp$ is the pseudo- 
scalar mass). For the y-axis we define the ratio of ro(jc)/ro(l), where x = 1 serves as a reference 



point to cancel the unknown overall scale. We have chosen to include QCDSF [ ]12[ ] and ETMC 



[|17[], because of the readily available data for ro and Mp$. The left panel shows the whole range of 
data points, whereas the right panel is a blow-up of the region close to the physical point. Assuming 
a physical value ro = 0.5 fm the physical value of x = r^Mp S is indicated by the vertical dashed 
line. As one can see the spread of the points is narrowing towards the physical point, where they 
start to agree within errors. 
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